The form of the elastic energy loss probability distribution in a static medium 
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We examine the probability distributions P(E,t) of the energy of a hard parton traveling in a 
partonic medium of constant density for a time t while undergoing elastic 2 2 pQCD interactions 
using a Monte-Carlo model. The form of these distributions is found to be non-Gaussian, confirming 
results by other groups with similarly detailed models and challenging the validity of the widely used 
diffusion approximation in elastic energy loss modeling. 
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I. INTRODUCTION 

One of the key observables in ultrarelativistic heavy ion 
experiments is the strong suppression of hadrons with 
high transverse momentum Pt in heavy-ion (A- A) col- 
lisions as compared with a p-p baseline. This has now 
been measured both at RHIC [1] and the LHC 0. This 
suppression is typically expressed in terms of the nuclear 
modification factor 



R AA (P T ,y,b) 



dNlJdPrdy 
T AA {h)daPP/dP T dy' 



(1) 



The transverse momentum dependence of R AA is notably 
more pronounced at the LHC than at RHIC, which is ex- 
pected on rather general grounds due to the harder par- 
tonic pt spectrum Q which can be understood as a filter 
through which the energy loss probability distribution 
P{AE) is observed. A more pronounced Pt dependence 
of R AA in turn implies a higher sensitivity to details of 
P(AE). Because of this increased sensitivity to the en- 
ergy loss probability distribution at the LHC, it is cru- 
cial for a model not only to reproduce the correct mean 
energy loss but also the correct fluctuations around the 
mean. 

The physics mechanism of the energy loss may be 
medium-induced radiation or elastic collisions with 
medium constituents. A sizeable elastic contribution to 
the total energy loss has been invoked by various groups 
[34To|. In such calculations, fluctuations of the elastic 
energy loss often are treated in a diffusion (Gaussian) 
approximation. In this approximation, a purely elastic 
picture reproduces the Py-dependence of Raa quite well 
and is in agreement with the LHC data [UlEI- However, 
this agreement is misleading, as in a calculation relaxing 
the diffusion assumption the form of R aa (Pt) resembles 
more the results of the radiative energy loss calculations 
and hence such a computation without diffusion approx- 
imation is in disagreement with the data [liij . In this 



paper, we illustrate that a large difference between the 
treatment in a phenomenological and a detailed Monte- 
Carlo (MC) model is the strongly non-Gaussian shape of 
P(AE). 

We carry out our study using a static QCD medium of 
constant temperature (and density). While unrealisti- 
cally simple, it nevertheless provides an useful setting for 
exploring the evolution of energy loss probability distri- 
butions, as well as providing a common ground on which 
to compare the results of different models. We first exam- 
ine how the shape of parton energy distribution depends 
on the assumptions made about the similarity of suc- 
cessive scatterings. We then compare the results of our 
model to those of Boltzmann Approach to MultiParton 
Scatterings (BAMPS) [HJ and results by Schenke et al. 

Hi- 



ll. THE MONTE CARLO SIMULATION 

Our energy loss model is based on the 2-to-2 scattering 
rates for a high-energy parton of a type i [T(j , 



where 



16tt 
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—2 / dE 2 f J (E 2 ,T)n iJ ^ kl (E 1 ,E 2 ,m 2 ), 

(2) 
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flij^.ki(Ei,E2,m 2 ) = ds[s(Tij-yki(s)]- (3) 

J2m 2 

The energy of the high-energy parton i is E\ and E 2 is 
the energy of the thermal quark (gluon) j with a Fermi- 
Dirac (Bose-Einstein) distribution function fj(E 2 ,T). 
The scattering cross section ay^jy (s) depends on the 
standard Mandelstam variable s. A regulator m = 
s m \/4ira s T is introduced for the infrared singularities 
appearing in the cross section. The strong coupling con- 
stant a s is kept fixed while s m is a parameter typically 
of the order of one. 

The propagation of the hard parton through the plasma 
is done in small time steps At. The Poisson probability 
for not colliding in this time interval is 
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P(No collisions in At) = e 



-r<Ai 



(4) 
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where I\ is the sum of all possible scattering rates for 
the hard parton of the type i. By keeping At small, we 
can safely assume that there will be at most one collision 
on this time interval. What type of scattering process 
happens is sampled according the ratios of the partial 
scattering rates ((2]) to the total scattering rate. 
In order to illustrate the effect of a diffusion approxi- 
mation, we use in this study two different methods for 
treating the momentum exchange between the hard par- 
ton and the medium particle. The standard one is to 
sample the Mandelstam-i from the % -distribution. The 
other method used here corresponds to the diffusion ap- 
proximation and requires us to calculate the average 
Mandelstam-i, 



(i) = 



s+m 2 dt 
s+m 2 al It 



(5) 



It is then assumed that all scattering processes happen 
with this average momentum exchange. Since (t) defined 
this way is not a constant but depends on Mandelstam- 
s, we additionally use the average values for the energy 
of a thermal parton (E2) = 3T and the collision angle 
(cos #12) ~ — I taken from [l^]. This leaves (s) dependent 
only on E\. In order to keep the approximation scheme 
simple, we also approximate the amplitude of the pro- 
cess qg — > qg with a pure f-channel: \M\ qg ^ qg oc s f 2 u . 
This approximation will decrease (t) slightly. This is of 
little consequence, as this study focuses on the variation 
in the form of the probability distributions and not in 
the magnitude of the average energy loss. With these 
approximations, (i) varies ~ 15% in the energy range 
Ex = 25 - 50 GeV at T = 300 MeV. 
In both schemes, after scattering the final state parton 
with highest energy is chosen as the new hard parton to 
be propagated further. The procedures outlined above 
are repeated for a predetermined time period t (corre- 
sponding to a length L), after which the simulation ends. 
Repeating the simulation several times produces a prob- 
ability distribution for the energy of the hard parton at 
the time t, P(E,t) ~ 1 dN(E ' t] 
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III. RESULTS 

We first examine the parton energy loss in a gluon 
plasma. Figure [1] shows the difference between the two 
scenarios described in the previous section for an initial 
50 GeV quark traveling in the medium for various times. 
Treating every collision as average event results in nearly 
Gaussian probability distributions for the parton energy, 
while in the scenario where we include the correct fluc- 
tuations by sampling the the thermal distributions and 
differential cross section the energy probabilities have a 
long tail towards lower energies. However, as we see in 
Figure [51 the dependence of the average energy loss on 
the in-medium pathlength L is linear in both cases. 



Next we select a gluon instead of a quark as initial par- 
ton and compare our results with the BAMPS [14j for a 
T = 400 MeV purely gluonic medium (i.e. we probe the 
channel process gg — > gg only) and using the same cou- 
pling constant value a s — 0.3. We find the qualitative 
behavior very similar, but quantitatively our model pro- 
duces stronger energy loss. This difference can be largely 
attributed to the differences in the cross section regular- 
ization, which according to Zapp et al. can produce a 
factor two difference in the cross section (TbT |. 
Finally we study a 10 GeV quark traveling through 
quark-gluon plasma with temperature T — 200 MeV, 
keeping the coupling strength a s = 0.3, and compare 
our results with those by Schenke et al. [15| • The results 
in this case agree both qualitatively and quantitatively. 
The main difference, i.e. the sudden drop in the 9.8-10.0 
GeV bin, is once again explained by the chosen regulator 
as we, unlike Schenke et al., ignore scattering events be- 
low the soft energy scale ~ ^/Air a s T and thus artificially 
suppress the number of particles with very small energy 
loss. Interestingly, our average energy loss is actually 
smaller than in Ref. UM in this scenario. 



IV. DISCUSSION 

We observe a strongly non-Gaussian shape of the energy 
(loss) probability distribution P(AE) of hard partons af- 
ter elastic collisions with a QCD medium of constant 
density. The non-Gaussian nature is mainly driven by 
fluctuations around the average energy transfer per col- 
lision, other effects like parton type changing conversion 
reactions with the medium (qq — > gg) which are often ne- 
glected in approximations to the elastic energy loss prob- 
lem do not play a role. 

Our results agree with observations previously made by 
Schenke et al. in Ref. fl5j . The emerging picture cast a 
serious doubt on the validity of the diffusion approxima- 
tion in the elastic energy loss modeling. 
The detailed functional form of P(AE) is not an aca- 
demic issue provided that the mean energy loss can be 
computed correctly, on the contrary it is directly linked 
to the Pt dependence of Raa at LHC kinematics. The 
implication of our results is that a large elastic contri- 
bution to the total energy loss (> 10%) is not indicated 
by the LHC data. Initially results using the diffusion 
approximation gave a reasonable agreement with the ob- 
served Pt dependence whereas radiative scenarios under- 
predicted the growth of Raa with Pt [ll| . This allowed 
for the possibility that a sufficiently large elastic contri- 
bution to the total energy loss could result in a scenario in 
agreement with the data. However, since a correct treat- 
ment of the fluctuations leads to the same underpredic- 
tion of the rise of Raa with Pt for elastic and radiative 
energy loss, this argument can no longer be made. 
We also find that the pathlength dependence of elastic 
energy loss remains linear (as expected) in a constant 
medium, no matter how fluctuations are treated, thus 
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FIG. 1: (Color online) a) Energy distribution of a 50 GeV quark after traveling in gluon plasma for various times, using average 
(t) (see text), b) Energy distribution of a 50 GeV quark after traveling in gluon plasma for various times, values of the t sampled 
from the cross section distribution. Dashed vertical lines illustrate the average energy values, while the solid line marks the 
initial energy bin. The temperature of the medium is T = 300 MeV, the strong coupling constant value is a a = 0.3 and the 
regulator parameter s m is unity. 



Path length dependence of energy loss 
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FIG. 2: (Color online) Comparison of the hard parton average 
energy pathlength dependence in average collisions scheme 
and the distribution sampling scheme. Solid circles corre- 
spond with the dashed lines shown in Fig. [T] a) and solid 
squares correspond with the dashed lines shown in Fig. Q]b). 
Temperature of the medium is T = 300 MeV, the strong cou- 
pling constant value is a 3 — 0.3 and the regulator parameter 
s m is unity. 



all arguments made previously in the context of a dif- 
fusion approximation (e.g. |17| |) based on pathlength to 
constrain the relative contribution of elastic energy loss 



FIG. 3: (Color online) Energy distribution of a 50 GeV gluon 
after various times, values of the t sampled from the cross 
section distribution. The vertical solid line marks the initial 
energy bin. The temperature of the gluonic medium is T = 
400 MeV, the strong coupling constant value is a s — 0.3 and 
the regulator mass parameter is s m = 1.0. Compare with 
upper panel of Fig. 4 in Ref. 



to be about 10% remain valid. Note also that dihadron 
correlations allow to constrain the magnitude of elastic 
energy loss from below and are consistent with a 10% con- 
tribution [18| . There is thus growing evidence that the 
relative contribution of elastic energy loss is considerably 
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FIG. 4: Energy distribution of a 10 GeV quark after traveling 
length L — 2 fm in a quark-gluon medium with a temperature 
T — 200 MeV, strong coupling constant value a a = 0.3 and 
mass parameter value s m — 1.0. The dashed vertical line 
illustrates the average energy and the solid line the initial 
energy bin. Compare with Fig. 3 in Ref. [l5| . 



smaller than expected by a straightforward computation 
based on an ideal quark-gluon gas picture. 
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FIG. 5: Energy distribution of a 10 GeV quark after traveling 
length L — 5 fm in a quark-gluon medium with a temperature 
T — 200 MeV, strong coupling constant value a s = 0.3 and 
mass parameter value s m = 1.0. The dashed vertical line 
illustrates the average energy and the solid line the initial 
energy bin. Compare with Fig. 4 in Ref. p^| . 
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